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ABSTRACT 

According to recent models, the accretion disk and black hole in active galactic nuclei are surrounded 
by a clumpy torus. We investigate the NIR flux variation of the torus in response to a UV flash 
for various geometries. Anisotropic illumination by the disk and the torus self-occultation contrast 
our study with earlier works. Both the waning effect of each clump and the torus self-occultation 
selectively reduce the emission from the region with a short delay. Therefore, the NIR delay depends 
on the viewing angle (where a more inclined angle leads to a longer delay) and the time response 
shows an asymmetric profile with a negative skewness, opposing to the results for optically thin tori. 
The range of the computed delay coincides with the observed one, suggesting that the viewing angle is 
primarily responsible for the scatter of the observed delay. We also propose that the red NIR-to-optical 
color of type-1.8/1.9 objects is caused by not only the dust extinction but also the intrinsically red 
color. Compared with the modest torus thickness, both a thick and a thin tori display the weaker NIR 
emission. A selection bias is thus expected such that NIR-selected AGNs tend to possess moderately 
thick tori. A thicker torus shows a narrower and more heavily skewed time profile, while a thin torus 
produces a rapid response. A super-Eddington accretion rate leads to a much weaker NIR emission 
due to the disk self-occultation and the disk truncation by the self-gravity. A long delay is expected 
from an optically thin and/or a largely misaligned torus. A very weak NIR emission, such as in hot- 
dust-poor active nuclei, can arise from a geometrically thin torus, a super-Eddington accretion rate 
or a slightly misaligned torus. 

Subject headings: accretion, accretion disks — dust, extinction — galaxies: active — galaxies: struc- 
ture — infrared: galaxies — infrared: ISM 
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1. INTRODUCTION 

Active galactic nuclei (AGNs) are powered by gas ac- 
cretion onto supermassive black holes (BHs) at the cen- 
ter of each galaxy. A variety of observations suggest 
that the accretion disk and the BH are surrounded by 
an optically and geometrically thick torus (Telesco et al. 
1984; Antonucci & Miller 1985; Miller & Goodrich 1990; 
Radovich et al. 1999). Since the torus potentially plays 
a role of a gas reservoir for the accretion disk, its nature, 
such as the structure, the size and the mass, has long 
been investigated (Pier & Krolik 1992, 1993; Fukue & 
Sanbuichi 1993; Granato & Danese 1994; Efstathiou & 
Rowan-Robinson 1995 Beckert k, Duschl 2004; Mor et al. 
2009). 

A large geometrical thickness of the torus revealed by 
various observations (Antonucci 1993; Pogge 1989; Wil- 
son & Tsvetanov 1994; Schmitt & Kinney 1996) indicates 
that numerous dusty clumps, rather than a smooth mix- 
ture of gas and dust, constitute the torus with a large 
clump-to-chmrp velocity dispersion ~ 100kms~ (Kro- 
lik & Begelman 1988; Wada & Norman 2002; Honig k 
Beckert 2007). Temperature of clumps is less than a 
critical temperature T su b ~ 1500 K above which dust 
grains are sublimated (Barvainis 1987). Infrared (IR) 
emission and absorption features provide unique oppor- 
tunities to probe the clumpy torus (Nenkova et al. 2002, 
2008; Dullemond & van Bemmel 2005; Honig et al. 2006; 
Geballe et al. 2006; Shirahata et al. 2007; Ibar & Lira 
2007; Schartmann et al. 2008; Deo et al. 2011). 

Clumps are heated by illumination from the central ac- 
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cretion disk, and closer clumps have higher temperature. 
Inner edge of the dusty torus is determined by the subli- 
mation process so that clumps' temperature equals T su b 
there, and radiates at Near-IR (NIR) as "3/im bump" 
(Rees et al. 1969; Neugebauer et al. 1979; Edelson & 
Malkan 1986; Kobayashi et al. 1993). Based on the en- 
ergy balance of the clump closest to the BH, Barvainis 
(1987) derived the innermost radius of the torus (dust 
sublimation radius, denoted as i? S ub,o m this study): 
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where Ltjv and a are UV luminosity and the size of dust 

grains, respectively. 

Indeed, NIR emission from type-1 AGNs lags behind 
optical variation by an order of a month (Clavel et al. 
1989; Glass 1992, 2004; Nelson 1996; Oknyanskij et al. 
1999; Minezaki et al. 2004; Suganuma et al. 2004). 
Moreover, the luminosity dependency of the time lag 
also coincides with the theoretical prediction as ex iyv 
(Suganuma et al. 2006; Gaskell et al. 2007). However, 
the NIR-to-optical time lag is systematically smaller 
than the lag predicted from Equation (1) by a factor 
of ~ 1/3 (Oknyanskij & Home 2001; Kishimoto et al. 
2007; Nenkova et al. 2008). To tackle with this conflict, 
Kawaguchi & Mori (2010, hereafter Paper I) pointed out 
that the illumination by an optically thick disk is in- 
evitably anisotropic, which is a fact missing in deriving 
Equation (1). There is a systematic difference between 
the inclination angle at which we observe the disk in type- 
1 AGNs and the angle at which an aligned torus observes 
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the disk. The effects of the anisotropic illumination nat- 
urally resolve the puzzle of the systematic deviation of a 
factor of ~l/3 (Paper I). 

In Paper I, we assumed the configuration appropriate 
for a typical type-1 AGN. In this study, we investigate 
the expected characteristics of the NIR emission for var- 
ious geometries of the disk, the torus and the observer. 
Anisotropic illumination by the disk and the effect of 
the torus self-occultation contrast our study with ear- 
lier works. The next section describes the calculational 
methods of our model. Then, properties of NIR emission 
from an aligned (Section 3) and a misaligned (Section 4) 
tori are presented. Finally, we make a summary of this 
study in Section 5. 

2. MODEL DESCRIPTION 

We calculate NIR reverberation/echo from the inner 
part of the torus in response to a flash of disk illumina- 
tion. The calculational method is essentially the same 
as Paper I, except for the incorporation of the torus self- 
occultation in this study. By considering the anisotropy 
of the disk illumination, we solved the inner structure 
of the torus, and explained why the observed time delay 
of NIR emission is systematically shorter than Equation 
(1). A large grain size and/or extinction between the 
torus and the disk are possibly alternative concepts to 
reduce the inner radius of the torus (Maiolino et al. 2001; 
Gaskell & Benker 2007; Gaskell et al. 2007; Kishimoto 
et al. 2007). 

While the geometry of the torus, the disk and the ob- 
server appropriate for a typical type-1 AGN was assumed 
in Paper I, we here investigate how the NIR response dif- 
fers with various possible geometries. A variety of type-1 
objects with a common i? S ub,o are compared. In other 
words, we compute for objects with the same isotropic- 
equivalent luminosity ("luminosity presuming isotropic 
emission") in optical/UV. 

Barvainis (1992) examined the NIR response due to 
the dust reverberation of an AGN torus. Since the inner 
radius of the torus is an input parameter there, the mean 
time delay of NIR emission behind the optical/UV flux 
variations is simply coupled with the assumed inner ra- 
dius. On the other hand, in Paper I and this study, both 
the inner radius and the time delay are output of calcu- 
lations. Moreover, he supposed that the whole torus, an 
ensemble of cube-shaped, optically thick clumps, is opti- 
cally thin. Since the inner part of an AGN clumpy torus 
is likely optically thick (Appendix), we take into account 
the torus self-occultation as well as the anisotropic emis- 
sion from spherical, optically thick clumps. Similarities 
and differences in the results for the time response be- 
tween our and his calculations are discussed in Sections 
2.4 and 3.1. 

In radiative transfer calculations of dusty clumpy tori 
(e.g., Honig et al. 2006; Nenkova et al. 2008), the 
isotropic illumination, mainly in optical and UV, is as- 
sumed for simplicity. We take into account the fact that 
the disk emission is inevitably anisotropic. Then, the 
broadband color between the optical/UV radiation from 
the disk and the NIR emission from the torus can be com- 
puted appropriately. In Sections 3.1 and 3.2, we compare 
our results with theirs. 

Below, our calculational method is summarised. 
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Fig. 1. — Radiation flux from the disk F(9) as a function of 
the polar angle 8, normalised to its pole-on value. Red solid line 
represents the net ^-dependency (projection plus limb darkening 
effects; Equation (2)), while black dashed line shows the former 
effect alone. 

2.1. Anisotropic Illumination of Disk 

Radiation flux (F) from a unit surface area of an op- 
tically thick disk toward a unit solid angle at the polar 
angle of 9 decreases with an increasing 9 as follows (Net- 
zer 1987): 

F{9) ex cos 6> (1 + 2 cos (9) (2) 

Here, the first term represents the change in the pro- 
jected surface area, while the latter represents the limb 
darkening effect for plasma, whose opacity is domi- 
nated by electron scattering over absorption (Sunyaev 
& Titarchuk 1985; Phillips & Meszaros 1986). Figure 1 
shows the ^-dependency of F(9), where the former ef- 
fect alone is also drawn for comparison. An accretion 
disk emits lesser radiation in the direction closer to its 
equatorial plane (i.e., larger 9; Laor & Netzer 1989; Sun 
& Malkan 1989; Hubeny et al. 2000). If the torus and 
the disk are aligned, the assumption of isotropic emission 
from accretion disks (e.g., Equation (1)) obviously over- 
estimates the radiation flux toward the torus, leading the 
overestimation of the inner radius of the torus. 

This effect works even if the disk is infinitesimally 
thin. As shown in Section 3.3, a nonzero thickness of 
the disk brings about another anisotropy of illumination 
flux, such that the torus is not illuminated below the 
disk height at 9 larger than a critical angle 9 max . Except 
in Section 3.3, we throughout adopt a thin disk with an 
aspect ratio of ~ 0.01, like the standard accretion disk 
model (Shakura & Sunyaev 1973). 

The effects of anisotropic emission and orientation have 
been discussed in the context of the Baldwin effect in 
the line fluxes of photo-ionised emission (Netzer 1985; 
Francis 1993; Bottorff et al. 1997). These effects upon 
the torus were examined for the first time in Paper I. 

2.2. Inner Structure of Torus 

The inner edge of the torus is determined so that the 
temperature of a clump (at the irradiated surface) equals 
T su b there. Since the radiation flux from the disk F varies 
with the polar angle 9, the sublimation radius of the 
torus R su b(9) is also a function of 9. Namely, R S ub(9) is 
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Fig. 2. — Schematic view of the inner structure of the torus, 
for a misaligned configuration (Section 4) with Storus ~ 20° and 
4'toTUB = 0°. Right cartoon indicates the definitions of # m i n (the 
torus thickness) and ma x (the disk thickness). 

the distance between the torus edge and the central BH 
for various 9. In contrast, we express the sublimation 
radius estimated under the isotropic emission assumption 
(Equation (1)) as i? S ub.o- The anisotropic illumination 
given in Equation (2) results in 
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Here, o b s is the polar angle toward the observer seen 
from the central accretion disk. Outside this radius, 
there are numerous clumps with their temperature be- 
low T su b- In the case of an isotropic emission from the 
disk, the torus edge is supposed to stand at a distance 
of i?sub.o- Although various grain sizes result in the sub- 
limation process occurring over a transition zone rather 
than a single distance (Nenkova et al. 2008), we employ 
a sharp boundary for simplicity. 

In Paper I, we found that (1) the torus inner edge is lo- 
cated closer to the central BH than suggested by previous 
estimations (Equation (1)) and that (2) the structure of 
the edge is concave/hollow. Moreover, (3) our result indi- 
cated that the innermost edge of the torus may connect 
with the outermost edge of the accretion disk continu- 
ously (e.g., Emmering et al. 1992; Elitzur & Shlosman 
2006). Figure 2 shows a schematic view of the torus in- 
ner region. If the torus is indeed a reservoir of gas for 
the disk, angular momenta of the infalling gas will align 
these axes. The misalignment between the torus and the 
disk is investigated in Section 4. 

2.3. Calculation of Transfer Function 

Current interferometric NIR observations (e.g., Swain 
et al. 2003; Kishimoto et al. 2009) and future Adap- 
tive Optics imaging with <~ 30 m telescopes cannot spa- 
tially resolve the innermost region (radius and shape) 
of the torus in nearby Seyfert galaxies. Thus, obser- 
vations of time variability will continue to be powerful 
tools to probe the innermost structure of the torus even 
in the coming decade. We calculate the time variation of 
NIR emission in response to a <5-function like variation 
of the optical/UV illumination [transfer function \&(£)]. 
Transfer functions contain various information of the re- 



emitting region, such as the shape and the emissivity pro- 
file etc (Blandford & McKee 1982; Netzer 1990; Gaskell 
et al. 2007). Time variation of the reprocessed radiation 
(NIR in this study) is a convolution of the illumination 
flux variation with ^(t). We calculate ^(t) and its cen- 
troid tdeiay, which corresponds to the observed time lag. 
Since the time profile is also one of the characteristics of 
this study compared to the earlier work (Sections 2.4 and 
2.5), we also calculate the width rms and the skewness 
s to describe the shape of $f(t): 
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«=/(*- idday) 3 *(t) dt/ *(t) dt J rms 3 (5) 

A negative (or positive) s means that the distribution has 
a left (or right) tail. The computed rms and s would 
be useful to interpret and predict the cross correlation 
functions between NIR and optical/UV light curves. A 
larger rms will correspond to a larger uncertainty in the 
measurements of the time delay. As to the detectability 
of NIR variations, a small rms and a high &(t) indicate 
the relative (i.e. in mag) and absolute (e.g., in ergs -1 ) 
variability, respectively. 

When the illumination flux varies, the inner edge of 
the torus shifts in principle (Laor 2004). Depending on 
whether dust grains in the clumps are sublimated or not, 
clumps belong to either the broad emission line region or 
the dusty torus (Netzer & Laor 1993; Suganuma et al. 
2006). However, it takes ~ lyear for the inner region 
of the torus to adjust to the varying illumination flux 
(Koshida et al. 2009; Pott et al. 2010). Thus, we regard 
that the inner structure of the dusty torus is steady in 
the timescale of NIR-to-optical time lag (~months). 

To calculate ^(t) for the clumpy torus, we consider 
the following items; (1) the optical path, (2) NIR emis- 
sivity of the torus inner region as a function of 9 and (3) 
anisotropic emission of each clump. In this work, (4) we 
include the effect of torus self-occultation (i.e., absorp- 
tion of NIR emission from a clump by other clumps on 
the line of sight). While the torus self-occultation is a mi- 
nor effect for a typical type-1 AGN, it plays a significant 
role for inclined viewing angles, thick tori and misaligned 
tori. Considering the self-occultation, we ignore the re- 
sponse from the aligned torus at 6 > \ (Appendix). 

First, (1) the optical path difference is written as 

-Rsub (0) [1 - {cos o bs cos + sin 9 ohs sin 9 cos cj>}] , (6) 

where (f> is the azimuthal angle and defined so that cf> = 
for the observer (Figure 2). The concave shape of the 
inner region of the torus reduces the optical path differ- 
ence. Clumps at slightly farther and those at slightly 
closer than i? S ub(0) will also emit NIR radiation, smear- 
ing out the resultant NIR response. Since this effect 
unlikely changes the ^(t) drastically, we consider only 
optical paths that hit the inner edge of the torus. 

Next, (2) for the emissivity as a function of 9, we as- 
sume that (2-1) the clump size increases and (2-2) the 
clump number density decreases when the clump-to-BH 
distance increases (e.g., Honig et al. 2006; Schartmann 
et al. 2008). Following Paper I, the emissivity of NIR 
flux per dfl is assumed to be proportional to i? S ub(0) 2 - 
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Fig. 3. — NIR response from 8 < ^ of the aligned torus, viewed 
from 9 ohs = 25° with t9 min = 45° and t9 max = 89°. Brown long- 
dashed line is computed by switching off both the waning effect 
and the torus self-occultation, making the left horn higher than 
the right one similar to the \l/(t) presented by Barvainis (1992). 
Purple short-dashed line then takes into account that each clump 
is optically thick, and is identical to the ^(t) presented in Paper 
I. Blue dotted line include both effects. Three solid lines show 
the contributions from each range. The portion at (a) </> ~ 
produces the red left bump, which is a faster (due to a shorter 
optical path difference) and weaker (due to self-occultation and 
a stronger waning effect) response than the opposite area at (c) 
(j> ~ 7T (green right horn). 

Third, (3) the anisotropy of the NIR emission from 
each clump is considered, since clumps are optically thick 
to NIR (and optical/UV) radiation (see Appendix for 
details). Namely, the question is how extent the illumi- 
nated surface of a clump is seen by the observer (Nenkova 
et al. 2002). Let us suppose that an observer looks at 
a clump with an angle £, where £ = means a face-on 
view of the illuminated surface. We adopt the following 
anisotropic coefficient for the waning effect, 
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This coefficient is chosen so as to reproduce the Monte 
Carlo calculations by Hdnig et al. (2006) for a single 
clump observed from three different £. 

Finally, (4) if the line of sight from a region to the 
observer passes through the torus, we omit the NIR flux 
from such a region. Radiation energy absorbed by the 
clumps on the way will be re-radiated at Mid-IR bands. 

The dotted line in Figure 3 presents the resultant 
transfer function of the aligned torus for t9 b s = 25°, cal- 
culated by summing up the NIR responses from different 
portions of the torus at (9, </>). The integration is done 
from t9 m i n to t9 max in the (9-direction and to 7r in the 
(^-direction. Here, the opening angle of the torus 8 min is 
assumed to be 45°, which is roughly consistent with vari- 
ous observational results (e.g., the ratio between types-1 
and -2 AGNs and/or opening angles of light cones in the 
narrow line region), and the maximum 9 of the torus 
#max is set to 89° (i.e., thin disk approximation). This 
parameter set is identical to the one adopted in Paper 
I, and is regarded as fiducial in this study. The mean 
delay ideiay, the width rms and the skewness s of ^(t) 



are 0.42 i? S ub,o/c, 0.13i? su b,o/c and —0.76, respectively. 
For typical Seyfert 1 galaxies and quasars (e.g., with UV 
luminosity of 10 43 ' 5 ergs _1 and 10 45 ' 5 ergs -1 ), the ex- 
pected inner radii -R S ub,o are ~0.1 and lpc (Equation 
1), while our resultant delays ideiay for the fiducial pa- 
rameter set are about a month and a year, respectively 
(see Figure 7). 

The long-dashed line is computed by switching off both 
the waning effect and the torus self-occultation, making 
the left horn higher than the right one similar to the time 
response of optically thin tori (Barvainis 1992). Then, 
the short-dashed line takes into account the waning ef- 
fect, and is identical to the \l/(£) presented in Paper I, 
with idciay, rms and s of 0.37i? su b,o/ c , 0.15i? su b,o/c and 
—0.35, respectively. The torus self-occultation results in 
a slight increase of ideiay and a more skewed profile. We 
will see the reason for these changes in the next subsec- 
tion. 

2.4. Response from various </> 

In order to clarify contributions from different (f>, we 
divide the torus inner edge into three regions equally. 
The three solid lines in Figure 3 show the ^(t) from each 
</> range. Among them, the left one is produced by the 
low <j) region (<j) = — -|), showing the rapid response 
at tdciay ~ 0.15i? su b,o/c (due to a short optical path 
difference) and a strong reduction of flux by the waning 
effect (large £) and the torus self-occultation. On the 
contrary, the right one with the longer delay at idciay ~ 
0.5i? su b,o/ c comes from the large </> region ((f) = ^- — 
7r), where clumps tend to direct their illuminated surface 
toward the observer, suffering from less waning effect. 
We see that the torus self-occultation affects no influence 
at <f) > $ with the fiducial parameter set. In other words, 
both the waning effect and the self-occultation selectively 
reduce the emission from the region with a short delay, 
making the right horn higher than the left one, opposing 
to the results for optically thin tori (Barvainis 1992). The 
fluences, j^(t)dt, from each <j> region are 1.3, 6.4 and 
10.4 from low to high 0, respectively. 

2.5. Response from various 9 

Later, we change the torus thickness (Section 3.2) and 
the disk thickness (Section 3.3). To understand how they 
will affect the time response of NIR emission, we draw 
\I/(£) from different 9 ranges (at 9 < 5) separately in 
Figure 4a. It turns out that most NIR flux arises from 
small 9 regions at 45° < 9 < 75° (i.e., the region with 
a high latitude from the equatorial plane), which show 
similar time delay. Because the emissivity is assumed 
in proportion to R su b(9) 2 , the fluence from the largest 
9 range (75°-89°) near the equatorial plane, where the 
torus inner edge is closest to the central BH, is small. 
The NIR fluences from each 9 range, from small to large 
9, is 8.4, 6.8 and 2.9, respectively. Therefore, little differ- 
ence in tdeiay for various t9 max and a shorter idciay for an 
extremely thin (0 m i n > 75°) torus are expected. Conse- 
quences of a thick torus are not drawn straightforwardly, 
because of the torus self-occultation. 

In case the torus is optically thin (cf. Appendix), we 
shortly mention the contribution from the far side (9 > 
4j) of the torus. Depending on the origin of optical/UV 
time variability of AGNs (X-ray reprocessing or change 
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Fig. 4. — (a) Transfer functions from different 9 ranges with 
the geometry same as Figure 3; 45°-60° (red solid line), 60°— 75° 
(green dotted line), and 75°-89° (black dot-dashed line). The NIR 
responses from 45° < 9 < 75° appear at similar time delay, (b) 
Transfer functions from the far side of the torus; 9 =91°— 105° 
(black dot-dashed line), 105°-120° (green dotted line), 120°-135° 
(red solid line), computed without the torus self-occultation. For 
comparison, ^(i) from 45°-89° (purple short-dashed line, as in 
Figure 3) is also drawn. If NIR emission from 9 > -^ penetrates 
the torus and arrives at the observer, tdclay will be longer. 

in the accretion rate etc.; e.g., Kawaguchi et al. 1998; 
Sakata et al. 2011), the link between the disk fluxes to 
9 > -| and to 8 < ^ will vary. The NIR time variations 
from 9 > ^ is controlled by the variations of the disk 
flux to 8 > 5, which we cannot observe. Thus, if the 
disk illumination to the two sides are random, the NIR 
flux variations from 9 > -| influence the measurements 
of the optical/UV-to-NIR lag as noise. On the other 
hand, if the time variations of disk illumination toward 
the two sides are similar, the time delay of NIR emission 
will become longer, as follows. To calculate the NIR 
response from 9 > f, we replace cos# in Equation (3) 
by | cos 6*|, and switched off the self-occultation effect. 
Figure 4b presents NIR responses from 9 > ^, which 
have long time delay due to their long light path. If 
no extinction affects the NIR emission from 9 > -| (as 
Barvainis 1992 assumed), we will obtain the net response 
from 45° < 9 < 135° with t dolay of 0.78i? sub ,o/c In other 
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Fig. 5. — Transfer functions for various viewing angles, with obs 
labelled near each line, from 0° (exactly pole-on geometry) to 44° 
(presumably corresponding to type-1.9 AGNs). Here, m i n and 
9 max are fixed at 45° and 89°, respectively. As 9 ^, s increases, (1) 
the centroid of the *(t) t dc i ay increases, (2) NIR flux [f^(t)dt] 
also increases, and (3) *&(£) becomes broad. (4) For nearly face-on 
geometry, \P(t) is peaky. 

words, a long ideiay may be a signature of an extremely 
low volume filling factor of clumps in the torus and/or a 
very thin torus (Appendix) . 

3. ALIGNED TORUS: VARIOUS DEPENDENCIES 

In this section, we present various dependencies of 
the NIR emission from the torus whose rotation axis is 
aligned to the disk axis. At the end of each subsection, 
we shortly summarise the obtained dependency by con- 
trasting with the result for the fiducial parameter set. 

3.1. Viewing Angle: # b s 

Figure 5 shows the transfer functions for various o bs, 
from an exactly pole-on geometry (0 o bs =0°) to inclined 
viewing angles. With a large (9 bs (« 40°-44°), the line 
of sight grazes the upper boundary of the torus, which 
would corresponds to the situation in type-1.8/1.9 AGNs. 
Here, 9 m [ n and 6> max are fixed at 45° and 89°, respectively. 
As o bs increases, we see that (1) the centroid of the re- 
sponse ideiay increases, (2) the NIR fluence J ^(t) dt also 
increases, and (3) the profile becomes broad. In addition, 
(4) ^(t) is quite peaky for a nearly face-on geometry. The 
first and second results are our new findings. Although 
the third and forth trends are already reported for op- 
tically thin tori by Barvainis (1992), we find here that 
both trends are also true for optically thick tori. 

First, (1) idoiay is drawn as a function of 9 b s in Figure 
6. For comparison, the result computed without the self- 
occultation is also shown. The torus veils selectively the 
region with a short delay (Figure 3), hence enlarges ideiay 
The self-occultation shows larger influences for more in- 
clined angles. On the right-hand side, we also draw a 
histogram of the observed delay in the unit of R sn b,o/c 
(based on Figure 7). The computed idciay ranges from 
0.27 to 0.74i? sub ,o/c (0.27 to 0.60i? sub ,o/c without the 
self-occultation), which covers the range of the observed 
time delay. In contrast, since Barvainis (1992) assumed 
an optically thin torus, ideiay was expected to have no 
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Fig. 6. — Centroid of \P(t) tdelay as a function of Q a \, a (left). 
Short-dashed line is computed without the torus self-occultation. 
While with this effect (solid line) , the computed tdelay ranges from 
0.27 to 0.74 -R su b,o/ c - It covers the observed range of tdelay shown 
in the right histogram based on Figure 7, where the typical error 
of tdelay is ~ 0.09 R S ub,o/ c - The viewing angle 6> t, s can be the pri- 
mary parameter to cause the observed scatter about the regression 
line in the tdelay - iuv diagram. 
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Fig. 7. — Relation between tdelay and the UV luminosity. Cor- 
responding scales are shown on the right axis. Observed data are 
taken from Suganuma et al. (2006) and references therein, as done 
by Kishimoto et al. (2007). Objects with multiple data are plot- 
ted by filled symbols; NGC5548 (blue squares), NGC4051 (red cir- 
cles), NGC7469 (purple inverted triangles) and NGC4151 (green 
diamonds). Three data points for NGC4151 are collected in 1969- 
2001, and are indicative of a long-term evolution of the viewing an- 
gle. Open symbols represent the objects with single data point for 
each: NGC3227 (triangle), Fairall 9 (diamond) , GQ Com (square), 
NGC3783 (circle) and Mrk744 (asterisk). Horizontal dashed lines 
show the ranges of the flux time variations. Dot-dashed line repre- 
sents -R su b o (Equation 1). Loci for various obs (from 0° to 44°) 
cover the observed scatter. 

(or quite weak) t9 b s -dependency. 

Since such a broad range is not achieved by the changes 
of min and # max as we will see later, we propose that the 
viewing angle is the key parameter responsible for the 
observed scatter about the regression line in the ideiay - 
Luv diagram (Oknyanskij & Home 2001; Suganuma et 
al. 2006). Conversely, the measurements of idciay are 
potentially useful to estimate the inclination angles. 

In order to compare our results with the observed data 
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Fig. 8. — Top: NIR fluence, f^(t)dt, as a function of t9 obs . As 
Figure 6, solid and short-dashed lines mean the results with and 
without the torus self-occultation. Dashed line indicates the disk 
flux toward = 0°, normalised so that it coincides with the NIR 
fluence at ohs = 0° . Middle: Width of the transfer functions rms 
(solid line) and the skewness s (blue dotted line). Inclined viewing 
angles lead to broader ^(t). Bottom: NIR fluence divided by rms 
(solid line) and peak value of the transfer functions (dashed line), 
which describe the peakiness of *&(t). At a nearly pole-on view 
with t9 t,s < 10°, *&(t) is peaky, implying the NlR variability is 
easily detected. 

more directly, the idciay v.s. £tjv diagram is drawn in 
Figure 7 (Suganuma et al. 2006 and references therein). 
Following Kishimoto et al. (2007), we estimate Luv by 
6 vL v {V). The uncertainty of the tdelay is ~ 0.09i? su b.o/c 
on average. Our loci for various # b s (with t9 m ; n and 
#max hxed at the fiducial values) well cover the observed 
scatter. A type-1.5 Seyfert galaxy NGC3227 (triangle) 
is located at an area with a small 6> b s (0°-30°), which 
would require a hysteresis effect (Koshida et al. 2009) 
or a thin torus for its relatively short tdelay Among the 
three points for NGC4151 (green diamonds), the lower 
two data (collected in 1969-1980 and 1990-1998) are con- 
sistent with a pole-on view, while the upper one in 2001 
indicates an inclined angle. A change of the viewing an- 
gle on the timescale of tens years, due to e.g., a precession 
of the disk, is indicated. 

Second, (2) the top panel of Figure 8 shows that the 
NIR fluence toward the observer, J \I/(t) dt, is insensitive 
to the viewing angle at 9 b s ^$ 25° and increases with # b s 
at more inclined angles (solid line). Comparing 9 ^ s = 0° 
with 44°, the fluence increases by a factor of 1.5. Due 
to the self-occultation, the NIR flux decreases by 19% at 
6» obs = 25° and by 26-28% at 6 ohs = 40°-44°. 

Regarding a continuous illumination from the disk as 
a series of flash, we can draw conclusions other than the 
time response of NIR emission. Namely, the obtained 
NIR fluence indicates the NIR flux toward the observer 
under a given steady optical/UV illumination. We recall 
that we compute for various objects with a common opti- 
cal flux toward the observer. Therefore, the y-axis (NIR 
fluence) directly indicates the NIR-to-optical color. In 
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other words, the NIR-to-optical color becomes red as 6* bs 
increases. We thus argue that the red IR-optical color 
observed for type-1.8/1.9 AGNs compared with type-1 
objects (Alonso-Herrero et al. 2003) arises not only from 
dust extinction but also from their intrinsically red color. 

The red color for large 9 bs objects can originate in the 
reduction of optical flux toward large 9 (Equation (2)) 
and/or in the change of NIR flux. Here, we try to distin- 
guish its origin. The dotted line shows the disk flux to- 
ward 9 — 0°, in proportion to [cos6> bs (1 + 2cos6> bs)]~\ 
normalised so that it coincides with the NIR fluence at 
#obs = 0°. If the NIR emission is isotropic, the NIR flu- 
ence in this diagram should change as the dotted line. In- 
stead, if the computed NIR fluence is smaller (or larger) 
than the dotted line, it means the NIR flux decreases (or 
increases) with 9. Our result indicates that the NIR flux 
decreases with 9 due to the torus self-occultation, but 
its ^-dependency is weaker than that of the disk flux, 
thereby presenting the red NIR-optical color for inclined 
angles. The reduction of the NIR flux toward a larger 
o bs is consistent with the radiative transfer calculations 
under an isotropic optical/UV illumination (Hbnig et 
al. 2006; Nenkova et al. 2008). Incorporation of the 
anisotropy of the optical/UV flux from the disk (Equa- 
tion (2)) enables us to predict the NIR-optical/UV color. 

Third, (3) the middle panel of Figure 8 describes the 
profile of \&(i). As 9 ^ s increases, the width of ^(t) 
rms also increases, meaning that ^(t) becomes broader. 
Both idciay and rms increase with # bs, and their ratio 
ttieizy/rms slightly rises toward a smaller viewing angle: 
a t&eUy/rms ratio increases twice between 9 b s = 40-44° 
and 0°. With a small o bs, the echo from various parts 
of the inner edge arrives at the observer at a similar de- 
lay, making the ^(t) quite narrow. For larger # bs, a 
variety of light pass difference arise between ~ and 
~ 7r, which causes the broader ^(t). We expect that the 
cross correlation function between NIR and optical/UV 
flux variations becomes broader for type-1. 5-1. 9 AGNs 
compared with a typical type-1 object. The skewness s 
is always negative (i.e., with a tail toward a shorter time 
delay), and the degree of asymmetry gets larger for a 
pole-on view. 

Finally, (4) we comment on the NIR variability ampli- 
tude (to be precise, the ratio of the NIR amplitude to 
the optical/UV one). A peaky $(£) will result in a large 
NIR amplitude, whereas a less peaky ^>(t) smears out 
the variability of reprocessed emission, producing a less 
NIR amplitude. In order to see this quantitatively, we 
draw the NIR fluence-to-rms ratio and the peak value of 
*(i) as a function of o bs in the bottom panel of Figure 
8. For nearly pole-on view with o b s < 10°, both quanti- 
ties rise. For a given optical/UV variability, such pole-on 
objects will show large NIR variability amplitudes. At 
$obs ^ 15°, the peakiness of \I/(t) is insensitive to the 
viewing angle. 

By contrast with the fiducial (9 b s of 25°, objects with a 
small viewing angle will exhibit a short time delay with 
a narrow and peaky response. On the other hand, a 
more inclined viewing angle leads to a longer delay with 
a broader profile and to an intrinsically redder NIR-to- 
optical color. The o b s -dependent delay contrasts our 
work with the earlier study for optically thin tori (Bar- 
vainis 1992). The computed range of tdeiay coincides with 
the observed one. The NIR response always shows an 
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Fig. 9. — Transfer functions for various torus thickness m i n . 
A smaller 6 m i n corresponds to a thicker torus. The torus self- 
occultation by the 9 < 60° region at </> ~ hides the left horn, 
which originates in the large 9 (60° < 9 < 80°) region at <j> ~ 
and is visible in green short-dashed line. 

asymmetry with a tail toward a shorter delay. 

3.2. Torus Thickness: 9 m i n 

In Paper I, we assumed that the semi-thickness of the 
torus (from the equatorial plane to the upper surface) is 
45°. However, luminous AGNs (quasars) seem to have 
thinner tori than Seyfert galaxies (Lawrence 1991; Ueda 
et al. 2003; La Franca et al. 2005; Arshakian 2005; Simp- 
son 2005; Maiolino et al. 2007; Hasinger 2008; Treister et 
al. 2008). Moreover, recent hard X-ray observations dis- 
covered type-2 AGNs with very thick tori (Levenson et 
al. 2002; Ueda et al. 2007; Eguchi et al. 2009; Noguchi 
et al. 2010). In this subsection, we show the expected 
characteristics of the NIR emission from type-1 AGNs 
with thick and thin tori. 

Figure 9 shows transfer functions for various 6 m { n , 
with # b s and # max fixed at 25° and 89°, respectively. 
When the torus thickens from 9 min = 45° to 26°, the 
self-occultation selectively veils the region with a short 
delay at <j) ~ 0, hence enlarges the delay and reduces 
the fluence, the width and s. In other words, the self- 
occultation by the region with <f> ~ and 9 < 60° hides 
the rapid response (left horn), which originates in the 
large 9 (60° < 9 < 80°) region and becomes visible when 
the torus gets thin (9 m m = 60°). On the other hand, the 
self-occultation unlikely influences the result for a thin 
torus. As shown in Section 2.5, a rapid response is ob- 
tained when the torus is extremely thin (# m i n > 75°), 
since R su b(9) is small at large 9. 

To see the # m i n -dependency more quantitatively, we 
draw idoiay as a function of 9 min in the top panel of 
Figure 10. Clearly, idciay becomes short for thin tori. 
Therefore, we expect that luminous quasars will show a 
relatively short time delay in the unit of i? S ub,o/c For 
instance, if the torus thickness is reduced for luminous 
objects, the loci in Figure 7 shown by the solid lines be- 
come bent (convex). The computed tdday ranges from 
0.16 to 0.47-R sub ,o/c for 6> min in 85°-26°. In contrast to 
the # bs-dependency, it seems difficult to explain the ob- 
served range (shown as the histogram in Figure 6) by a 
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Fig. 10. — Top: Time delay (in the unit of ii au b o/ c ) as a function 
of torus thickness £? m i n . The meanings of different lines are the 
same as Figures 6 and 8. A thick torus veils the region with a short 
light path difference at <j> ~ 0, and hence enlarges idelay an d reduces 
rms at 6? m i n < 55°. As the torus becomes thinner, tdelay becomes 
shorter. Middle and Bottom: The same as Figure 8, but for the 
0min-dependency here. As the torus becomes thicker, \&(i) shows a 
longer delay with a narrower and more heavily skewed profile due 
to the torus self-occultation. A less peaky \P(t) is expected for a 
thin torus, making the detectability of NIR variability difficult. 

change of 9 m [ n alone. 

As the torus becomes thicker, ^(t) shows a longer 
delay with a narrower and more heavily skewed profile 
at 9 m [ n < 55° due to the torus self-occultation (mid- 
dle panel). This trend opposes to the results for various 
viewing angles, where tdelay and rms are positively cor- 
related each other. Although a low ^>(t) is expected for a 
thin torus (bottom panel) , the small rms implies a large 
relative variability in the NIR emission. 

Next, Figure 9 also shows that the NIR fluence is a 
strong function of 9 min . Qualitatively, it is trivial, since 
various 9 min mean various solid angles of the torus sub- 
tended at the central BH, 51 t0 rus- Here, 
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However, as shown in numerical results by Nenkova et 
al. (2008), the NIR flux is not exactly in proportion to 
fitorus- [Therefore, we need a caution when relating the 
observed NIR-to-UV luminosity ratio with floras (e.g., 
Mor & Trakhtenbrot 2011; see also Section 3.3).] More- 
over, we let the inner radius of the torus R su h{9) vary 
with 9. Thus, it is not obvious how the NIR fluence of 
our torus model varies with ri t0 rus- 

Figure 11 shows that the NIR fluence decreases dras- 
tically as 9 m in increases. A thin torus locates its inner 
radius at a short distance from the BH (Equation (3) 
and Figure 2), where the size of clumps is small and the 
NIR emissivity is low (Section 2.3). For thin tori with 
$min ^ 55°, we find that the NIR fluence is roughly in 
proportion to fit r US - Thus, luminous quasars are ex- 
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Fig. 11. — NIR fluence, f^(t)dt, to the observer at 6> obs = 25° 
for various torus thickness 6 min , as a function of the solid angle of 
the torus seen from the central BH f2torus. Solid and short-dashed 
lines have the same meanings as Figures 6, 8 and 10 top panels. For 
the former, circles are also plotted in 5° steps of m i n , and labelled 
in 10° steps. As the torus becomes thicker with m ; n < 40°, the 
NIR flux starts to decreases owing to the torus self-occultation. 
For a thin torus with t9 m i n > 55°, the NIR fluence is in proportion 
to £^J j US - Green dotted line with triangles is based on the 2/im 
flux densities calculated by Nenkova et al. (2008), for the torus 
semi-thickness parameter cr of 15°, 30°, 45° and 60° (from the left 
to the right). A modest thickness of the torus leads to the strongest 
NIR emission. 

pected to show weak NIR emission (blue NIR-to-optical 
color), which is consistent with the observed trend of the 
decreasing NIR-to-optical flux ratio with an increasing 
optical luminosity (Maiolino et al. 2007; Treister et al. 
2008; Jiang et al. 2010; Mor et al. 2011). AGNs with 
very weak NIR emission (such as " hot-dust-poor" AGNs 
named by Hao et al. 2010) may indicate that their tori 
are very thin. For comparison, we also draw the flux 
densities at 2/im from clumpy tori seen from 9 \> s = 25°, 
based on a radiative transfer computation (Figure 8 of 
Nenkova et al. 2008, corrected following its Erratum by 
Nenkova et al. 2010, and scaled to match with our result 
at t9 min ~ 45°). They present the results for four different 
torus semi-thickness cr, with a gaussian clump distribu- 
tion. For the smooth boundaries of the torus in the 9 
direction, the solid angle is computed as (see Apendix): 
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which is larger than cos (f — cr). The steeper decline of 
the NIR flux in our result is likely due to the 9 depen- 
dency of i? S ubi($)> hi contrast to the constant sublimation 
radius (resulted from the presumed istropic illumination) 
in their computations. 

Next, we move on to the result for a thick torus. 
The self-occultation reduces the NIR flux by 55-67% for 
9 m in = 30°-26°, whereas it results in a reduction by 19% 
for the fiducial 9 m [ n of 45°. Due to the self-occultation, 
too thick tori with t9 min < 40° exhibit lesser NIR flux 
(i.e., intrinsically bluer NIR-optical color) than a torus 
with a moderate thickness with 6> m j n ~ 40°-45°. The 
flux reduction of a thick torus compared with a mod- 
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assumed. 



As mentioned in Section 2.5, idelay seems insensitive to 



estly thick torus is again consistent with the result by 
Nenkovaet al. (2008). 

Bringing together the two behaviours, both a thick and 
a thin tori show the weak NIR emission. Namely, a mod- 
est thickness of the torus leads to the strongest NIR emis- 
sion. Therefore, a selection bias will arise such that NIR- 
selected AGNs tend to possess moderately thick tori. 

Compared with the modest torus thickness (# m i n = 
45°), both a thick and a thin tori display the weaker NIR 
emission (i.e., a bluer NIR-optical color). A thick torus 
shows a slightly delayed, narrow and largely skewed NIR 
response. On the other hand, as the torus gets thinner, 
the NIR response becomes more rapid, narrower, closer 
to time-symmetric and low. 

3.3. Disk Thickness: # max 

When the accretion rate exceeds the Eddington rate 
(« 16I/Edd/c 2 ), an optically thick advection-dominated 
accretion flow (a slim disk) appears (Abramowicz et al. 
1988). Since super-Eddington disks are geometrically 
thick (Abramowicz et al. 1988; Madau 1988), they can- 
not illuminate the directions near their equatorial plane 
by the disk self-occultation (Fukue 2000). 

As discussed in Paper I, some AGNs with presum- 
ably super-Eddington accretion rates show the weak 
NIR emission (Ark564, TonS180, J0005 and J0303; 
Rodriguez- Ardila & Mazzalay 2006; Kawaguchi et al. 
2004; Jiang et al. 2010; see, however, Hao et al. 2010). 
Small # max due to the self-occultation by a geometri- 
cally thick disk can be a reason for the weakness. More- 
over, the observed data do not support the concept of 
Eddington- limited accretion (Collin & Kawaguchi 2004). 
Thus, a strong anisotropy of the disk emission, such as 
Equation (2) and the disk self-occultation, is required 
to allow gas infall to super-Eddington accreting sources. 
In this subsection, we examine the influences of the disk 
thickness and the accretion rate upon the NIR emission 
in more detail. 

In principle, the disk thickness is a function of the 
distance from the central BH. Both the illumination 
spectrum from an AGN disk and the absorption effi- 
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Fig. 13. — Same as Figures 8 and 10, but for the 6 m ax- 
dependency. Corresponding gas accretion rates onto the BH in the 
unit of -LEdd/c 2 are labelled in the top panel for tdelay Second: 
Larger accretion rates make the disk thicker (i.e., smaller <9 max ) and 
the shade of the disk itself larger (i.e., smaller ntorus). Thus, the 
NIR fiuence becomes small as the accretion rate increases. Third: 
rms and s as well as tdelay are i nsens itive to max . Bottom: With 
high accretion rates, the NIR variability will be relatively small. 



ciency of dust have their peaks at Far-UV (e.g., Laor 
& Draine 1993). Therefore, we deduce the disk thick- 
ness at the Far-UV emitting region. Based on the work 
by Kawaguchi (2003; his Figure 5 for a 10 65 M BH), 
the semi-thickness of the disk (=90°— 6 max ) at the re- 
gion with the temperature of (4-5)10 4 K are 1°, 4°, 
17° and 39°, for the accretion rates of 1, 10, 100 and 
1000LEdd/c 2 , respectively. Comparing the first and the 
last cases, the solid angle of the torus illuminated by the 
disk (Equation 8) differs by a factor of ~ 9 (a smaller 
^toms with a higher accretion rate). 

Transfer functions for four accretion rates are pre- 
sented in Figure 12. As we see in Section 2.5, tdelay is 
insensitive to # max (top panel of Figure 13). Therefore, 
the observed broad range of tdciay is not reproduced by 
a change of the accretion rate. Larger accretion rates 
make the disk thicker and the shade of the disk itself 
larger (i.e., smaller ^toru S ; Equation (8)). Thus, the NIR 
fiuence becomes small as the accretion rate increases, as 
shown in the second panel. By changing the accretion 
rate from 1 to 1000LEdd/c 2 , the fiuence becomes ~ 1/5. 
This is consistent with the weakness of the X-ray emis- 
sion line from neutral iron, which potentially originates 
in the illuminated torus, in a narrow-line quasar (Taka- 
hashi et al. 2010; cf. Page et al. 2004). Similar to tdelay, 
rms and s are insensitive to # max (third panel). On the 
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contrary, the height of &(t) (bottom panel) is affected 
by the accretion rates in the sense that super-Eddington 
sources will show a low absolute variation. 

Other than the reduction of 9 max and fitoms, super- 
Eddington accretion rates cause another influence upon 
the NIR flux. Sub-Eddington accretion disks suffer little 
from the disk self-gravity, hence extend far away from the 
central BH, radiating across UV, optical and NIR bands 
(Tomita et al. 2006; Kishimoto et al. 2008). When the 
accretion rate becomes super-Eddington, the disk self- 
gravity starts to govern the disk and truncate the outer 
part of the disk. Due to the truncation, super-Eddington 
disks do not radiate at NIR (Kawaguchi et al. 2004). 
To sum up, both small max of the torus and the small 
outer radius of the disk, caused by the high accretion 
rate, provide less NIR emission. 

In contrast to the result for a thin disk with a sub- 
Eddington accretion rate of ^Edd/c 2 , a super-Eddington 
accretion rate leads to a much weaker NIR emission (with 
a help of the disk self-gravity) and to a low time response. 

Here, we summarise the three dependencies examined 
above. Only the variation of the viewing angles repro- 
duces the observed range of the time delay. Therefore, 
we propose that the viewing angle is the key parameter 
responsible for the observed scatter in the ideiay-^uv di- 
agram. A weak NIR emission (such as in hot-dust-poor 
AGNs in the weakest cases) indicates either a thin torus, 
a thick torus or the super-Eddington accretion. On the 
other hand, a stronger NIR flux (a redder NIR-optical 
color) than the fiducial one is obtained only by a large 
viewing angle for a modestly thick torus. 

4. MISALIGNED TORUS 

So far, we have assumed that the rotation axis of the 
torus is aligned to that of the disk. Since a specific an- 
gular momentum is likely larger in the torus than in the 
disk, the assumption sounds plausible. Indeed, we have 
examples, where the rotation axes of maser disks are 
aligned to the jet directions (e.g., Miyoshi et al. 1995; 
Meisenheimer et al. 2007; Mamyoda et al. 2009). How- 
ever, there are also a number of examples showing the 
misalignment between the maser disks and the jets (e.g., 
Yamauchi et al. 2004; Raban et al. 2009). Random ori- 
entation between jets and galactic axes is also reported 
(Clarke et al. 1998; Nagar & Wilson 1999; Schmitt et al. 
2001). 

In this section, we investigate the consequences of the 
misalignment between the torus and the disk axes. Now, 
we introduce #toms and 0toms to specify the rotation axis 
of the torus relative to that of the disk and the observer 
(0 = 0). Figure 2 shows the geometry with #torus ~ 20° 
and chorus = 0°. 

Three #toms are examined; 30°, 60° and 90°. In the 
first two cases, we adopt o b s of 25°. Since the direction 
to # bs of 25° is obscured by the torus if #torus = 90°, 
#obs of 65° is chosen in the last case. With t orus = 60° 
and 90°, the inner edge of the torus at 9 > ^ becomes 
visible. Thus, for these two #torus, we integrate 9 not only 
from to # m ax but also from 7r — # max to 7r. Here, the 
semi-thickness of the torus and the disk are fixed at 45° 
and 1° (i.e., 9 m [ n = 45° and 9 max = 89°), respectively. 

Figure 14 shows the transfer functions (integrating cf> 
from to 27r), and Table 1 summarises the results of a 
series of calculations. For reference, we also show the 
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Fig. 14. — Transfer functions for the misaligned tori with various 
^torus and </»torus- For reference, open blue circles indicate the 
response of the aligned torus seen from f3 b s = 25° (identical to 
the twice enlarged version of the blue dotted line in Figure 3). (a) 
0torus = 30° and (b) 6> t orus = 60° with f? obs = 25°. (c) 6» torus = 90° 
with f? obs = 65°. In (b) and (c), the NIR emission from 8 > 7r/2, 
which produces the NIR response at a long time delay, is visible. 

response of the aligned torus from < <fi < 2ir, which is 
identical to the twice enlarged version of the dotted line 
in Figure 3 (obtained by < < 7r integration). For 
#torus of 90°, the disk illumination flux toward the torus 
is larger than the flux toward the observer, meaning that 
the inner radius of the torus is large. Blank fields (and 
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TABLE 1 
Results for Misaligned Tori. 
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0.28 

0.003 
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44 
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Note. - Blank fields (and </>torus > 120°) mean that the 
direction to the observer (at <j> = 0) from the central BH is obscured 
by the torus with such parameter sets. 

a Hereafter, the time delay is quoted in this unit. 

b NIR flucncc for the fiducial parameter set is calculated by in- 
tegrating 9 from to 27T, which has twice the fluence shown in 
Sections 2 and 3. 

'/'torus > 120°) in the Table mean that the direction to 
the observer from the central BH is obscured by the torus 
with such parameter sets. In other words, a large </>torus 
excludes type-1 AGNs, and is biased to type-2 AGNs. 

A variety of tdciay is achieved from 0.25 to 
0.98.R su b,o/c, which is wide enough to match up with 
the observed range (Figure 6). A short tdeiay is obtained 
for 0torus of 30°, and is associated with a small NIR flu- 
ence. The geometry with 6>torus = 30° and </>torus = 0° 
is similar to a pole-on view of aligned tori (Section 3.1), 
producing a very narrow NIR response at a short time 
delay. On the other hand, a large tdeiay is realized for 
#torus of 90°, accompanied with a fluence similar to the 
one for the fiducial parameter. 

As the torus becomes inclined relative to the disk, the 
torus hides not only the region with a short delay but 
also various area of the inner edge of the torus, making 
the NIR fluence small and the time response complicated. 
Contributions from 9 > ^ , which are observable for #torus 
of 60° and 90°, appear at a longer time delay. Thus, they 
tend to increase the time delay, the NIR fluence and the 
skewness (s > 0). With 6> t0 rus = 90°, the NIR emis- 
sion from 9 > ^ is large enough to compensate the flux 
reduction due to the torus self-occultation, exhibiting a 
comparable fluence to the reference result. A long idciay 
(~ -RsubjO/c), if observed in the future, may be a signa- 
ture of a large # t orus (^90°; i.e., heavily misaligned). 

In case the torus self-occultation does not work due 
to an extremely small volume filling factor of clumps etc 
(Appendix), an aligned torus shows the NIR response 
as described in Section 2.5 with idciay = 0.78i? SUD ,o/ c - 
Summing up all the four lines in Figure 4b and multiply 
the fluence by two (to convert the < <f> < it integration 
to the 0-2-7T integration), we obtain the fluence of 140. In 
contrast, we find that misaligned tori without the self- 
occultation effect show a yet longer delay with a huge 
NIR fluence. For instance, for 6> t0 rus = 60° with 9 ^ s = 



25° and for 6> t0 rus = 90° with 9 b s — 65°, idciay = 1.2 and 
1.9 i?sub,o/c with the fluence of 250 and 960, respectively. 
Therefore, an extremely long (> -R su b,o/ c ) delay with 
a big NIR flux would mean a misaligned optically thin 
torus. 

5. SUMMARY 

According to recent models, the accretion disk and the 
BH in AGNs are surrounded by a clumpy torus, with its 
inner radius governed by the dust sublimation process. 
Regarding the inner radius of the torus, there was a sys- 
tematic deviation between the observational results and 
the theory. In Paper I, we showed that the anisotropy of 
the disk emission resolves this conflict for a typical type- 
1 AGN. We found that the anisotropy makes the torus 
inner region closer to the central BH and concave. Fur- 
thermore, the innermost edge of the torus may connect 
with the outermost edge of the accretion disk continu- 
ously. 

In this study, we have calculated the NIR flux variation 
of the torus in response to a UV flash for various geome- 
tries of the disk, the torus and the observer. Anisotropic 
illumination by the disk and the effect of the torus self- 
occultation contrast our study with earlier works. We 
have found that both the waning effect of each clump and 
the torus self-occultation selectively reduce the emission 
from the region with a short delay. Thus, the resultant 
NIR time response shows a # ODS -dependent delay and an 
asymmetric profile with a negative skewness, opposing to 
the results for optically thin tori (Barvainis 1992). 

By contrast with the fiducial viewing angle of 25°, a 
small viewing angle results in a short time delay with 
a narrow and peaky response. On the other hand, a 
more inclined viewing angle leads to a longer delay with a 
broader profile and to a redder NIR-to-optical color. We 
propose that the red NIR-optical color of type-1. 8/1. 9 
objects is caused by not only the dust extinction but 
also intrinsically red color. The computed range of idciay 
coincides with the observed one. 

Compared with the modest torus thickness of 45°, both 
a thick and a thin tori display the weaker NIR emission, 
consistent with the work by Nenkova et al. (2008). In 
other words, a modest thickness of the torus leads to the 
strongest NIR emission. A selection bias is thus expected 
such that NIR-selected AGNs tend to possess moderately 
thick tori. For thin tori, we have found that the NIR flu- 
ence is in proportion to f2^ us . As the torus becomes 
thicker, the NIR response shows a slightly longer de- 
lay with a narrower and more heavily skewed profile due 
to the torus self-occultation. This trend opposes to the 
viewing angle dependency, where the delay and the width 
are positively correlated each other. On the contrary, as 
the torus gets thinner, the NIR response becomes more 
rapid, narrower, closer to time-symmetric and low. 

In contrast to a thin disk with a sub-Eddington ac- 
cretion rate, a supcr-Eddington accretion rate leads to a 
much weaker (< 1/5) NIR emission due to the disk self- 
occultation and the disk truncation by the self-gravity, 
and to a low time response. 

Among the three dependencies examined for aligned 
tori, only the variation of the viewing angles reproduces 
the observed range of the delay. Therefore, we propose 
that the viewing angle is the key parameter responsible 
for the observed scatter about the regression line in the 
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idoiay v.s. Luv diagram. Conversely, the measurements 
of ideiay are potentially useful to estimate the inclination 
angles. 

We have also investigated the consequences of the mis- 
alignment between the torus and the disk axes. A vari- 
ety of ideiay is achieved, which is wide enough to cover 
the observed range. A short delay is obtained for a small 
misalignment and is associated with a small NIR flucnce, 
while a a long delay is obtained for a largely misaligned 
torus with an usual fluence. This trend contrasts with 
the viewing angle dependency for aligned tori, where a 
short delay is associated with a normal NIR fluence while 
a long delay means a large NIR fluence. For highly mis- 
aligned cases, contributions from 9 > ^ increase the time 
delay (up to ~ R su b.o/c), the NIR fluence and the skew- 
ness (s > 0). 

In case the torus is optically thin (with an inefficient 
self-occultation) , the time delay of the NIR emission from 
an aligned torus becomes longer. Moreover, misaligned 
optically thin tori show a yet longer delay (> i? S ub,o/c) 
with a huge NIR flux. 

From an observational point of view, these numerical 
results are summarised as follows. If the observed time 
delay of the NIR emission is short, it will mean either 



a small viewing angle, a geometrically thin torus, or a 
slightly misaligned torus. If the delay is long, on the 
other hand, it indicates either an inclined viewing angle, 
an aligned optically thin torus, or a largely misaligned 
torus. An extremely long delay (> i? su b,o/c) would mean 
a misaligned optically thin torus. As to the NIR flux, a 
blue NIR-to-optical color (i.e., a weak NIR emission such 
as in hot-dust-poor AGNs in the weakest cases) indicates 
either a geometrically thin torus, a geometrically thick 
torus, a super-Eddington accretion, or a slight misalign- 
ment between the torus and the disk. On the other hand, 
a red NIR-optical color (a large NIR emission) means a 
large viewing angle with a modest geometrical thickness 
of the torus or a largely misaligned optically thin torus. 
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APPENDIX 
OPTICAL THICKNESS OF CLUMPY TORI 

We briefly summarise the basic properties of clumps in the torus. The main aim here is to deduce the optical 
thickness of the clumpy torus. If it is optically thin at NIR (as Barvainis 1992 assumed), we need to consider the 
NIR emission from all the clumps including those at the far side of the torus which is illuminated by the back side 
(relative to the observer) of the accretion disk. In contrast, if it is optically thick, the torus self-occultation should be 
incorporated appropriately. 

Each clump located at a distance of r from the central BH is characterised by the radius of i? c iump and the mass of 
M c i ump . Clumps that survive against its own thermal pressure and tidal shearing by the BH must be heavier than the 
Jeans mass and/or more compact than the tidal (Hill or Roche) radius (e.g., Vollmer et al. 2004). Marginally stable 
clumps, at the boundary of these criteria, have the following various quantities (Hdnig & Beckert 2007), 
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Here, p c iump, Ah, clump, c s and m p are the mean density in the clump, the column density of the clump, the sound 
speed in the clump and the proton mass, respectively. The normalization in p c iump (corresponding to 310 8 cm~ 3 ) is 
consistent with the observed lower limit for the mean number density (>10 7 cm -3 ; Geballe et al. 2006; Shirahata et 
al. 2007). Incidentally, clumps in the broad line region (at r ~ 10 16 cm) have their sizes around 10 13 cm (Risaliti et al. 
2009; Maiolino et al. 2010), consistent with the extrapolation of Equation (Al) to smaller r. However, their density 
(~ 10 n cm~ 3 ) and iVn, clump [~ (2 — 9) 10 23 cm~ 2 ] are much less than those expected from Equations (A3) and (A4). 



The column density above means Compton-thick (~ 5 a, 



Thomson)' anc ^ corresponds to optical depths at V- and 
-ftT-bands of 1400 and 160, respectively, using the conventional extinction law (Savage & Mathis 1979; Cardelli et 
al. 1989). In radiative transfer calculations of AGN tori, clumps with optical depth at y-band of 30-100 are often 
adopted (Ncnkova et al. 2008; Honig et al. 2008; Deo et al. 2011). Krolik & Bcgclman (1988) estimates A^h, clump to 
be ~710 23 cm~ 2 (^one tenth of Equation (A4)). Even with such a reduction of the clump opacity by jn—h;, each 
clump is still opaque to NIR photons. 

Therefore, optical thickness of the torus r t0 rus is simply related to the probability that the incoming rays hit clumps. 
We treat the two directions separately: one in the vertical direction (parallel to the rotation axis) with a suffix of 
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Fig. 15. — Optical thickness of the clumpy torus in the vertical direction as a function of the distance from the central BH r. Red solid 
and black dotted lines are for Mbh of 1O 7 M0, with the clump radius assumed to be one tenth of that in Equation (Al) in the latter. Blue 
dashed line represents the case of 10 9 Mq. For all the cases, the inner part of the clumpy torus is optically thick. 

||, and another along the equatorial plane with a suffix of _L. The key parameter here is the volume filling factor of 
clumps in the torus /, for which we assume 0.03 (Vollmer et al. 2004). By writing the number density of clumps in 
the torus by n c , 
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where we assume the thickness of the torus H to be ~ r. Now, T torus means the mean number of clumps along the 
ray, and a fraction e _Ttorua of the incoming rays pass through the torus without encountering any clumps (Natta & 
Panagia 1984; Nenkova et al. 2002). 
If we adopt the clump size in Equation (Al), then 
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Assuming that clumps are mainly heated by the direct illumination from the central accretion disk (i.e. ignoring 
the irradiation from nearby clumps), the clump temperature will be proportional to r -°- 5 r, - 25 ^ with L being the 
illumination luminosity (cf. Nenkova et al. 2008). Therefore, the size of torus emission is larger for larger L at longer 
A, in the form ex X 2 L 05 (Tristram et al. 2009). Then, we assume for c s as follows, 
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with LEdd being the Eddington luminosity. In summary, 
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Figure 15 shows the optical thickness of the clumpy torus as a function of r. Throughout, we fix Co, L/L^dd and / so 
that these in parentheses in Equation (A10) equal unity. For instance, T torus n is about five at r of O.lpc. 
The radius where the torus is opaque r opa quR is achieved by setting Tt nrns h = 1: 
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The inner part within this radius of the torus is likely optically thick. In this study, we therefore consider basically 
emission from the near side of the torus, as we did in Paper I. This assumption can be tested in principle via the 
profile of broad emission lines and its time variation (e.g., Peterson 2001). If some indications of the emission from 
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the far side of the torus is observed, the volume filling factor seems extremely low (<C 0.03/5) or the torus is very thin 
(F/r<l/5). 

Similarly, we deduce the optical thickness perpendicular to the rotation axis Ttorus._i_j which is equivalent to Afo (the 
average number of clumps along radial equatorial rays) denoted by Nenkova et al. (2008). They showed that Ao is 
likely between 5 and 15, consistent with our estimation below. Adopting the inner radius of i? su b,o an d the clump size 
in Equation (Al), 
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Here, we assume L = 2.5 L\jy. Since this optical thickness is also larger than unity, we restrict ourselves to non-obscured 
objects (i.e., the geometry where the line of sight to the central BH is not blocked by the torus). 
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